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We address the properties of surface-wave sohtons on the interface between a semi-infinite homo- 
geneous linear medium and a semi-infinite homogeneous nonlinear nonlocal medium. The stability, 
energy fiow and FWHM of the surface wave solitons can be affected by the degree of nonlocality of 
, the nonlinear medium. We find that the refractive index difference affects the power distribution 

04 . of the surface solitons in two media. We show that the different boundary values at the interface 

< I ' can lead to the different peak position of the surface solitons, but it can not infiuence the solitons 

' stability with a certain degree of nonlocality. 
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PACS numbers: 42.65-k, 42.65.Tg, 42.70.Df 

Surface waves propagating along the interface between a homogeneous linear medium and a homogeneous nonlinear 
^ medium display many interesting properties, which have no analogues in homogeneous media. Decades years ago, 
. ^ _ such surface waves had been studied in the local nonlinear optical case [H-Q . [1] and [6] show that there is no stable 
^ • surface wave when the zero field refractive index of the nonlinear medium is larger than the refractive index of the 
Q ' linear medium, on the contrary, surface wave is stable. 

^ In the nonlocal nonlinear optical domain, such surface waves were analyzed at the interfaces of diffusive Kerr-type 

materials f7"9] . Recently, surface- wave solitons were observed at the interface between a dielectric medium(air) and a 
nonlocal nonlinear medium(lead glasses) They found that these solitons are always attracted toward the surface, 
^>-^ and unlike their Kerr-like counterparts, they do not exhibit a power threshold. Two-dimensional surface solitons 
] featuring topologically complex shapes, including vortices and dipoles with nodal lines perpendicular to the interface 
Oh. of nonlocal thermal media were studied in [11]. Defocusing thermal materials can also support surface waves under 
appropriate conditions jH, [l3|- Multiploe solitons localized at a thermally insulating interface are addressed in [14]. 
^-H . However, to our knowledge, the variation of such surface- wave solitons due to the change of the degree of nonlocality 
^ ■ or the boundary value at the interface of the semi-infinite nonlocal nonlinear media and the semi-infinite linear media 
, were not studied to this day. In this Letter, we reveal that the degree of nonlocality can affect the stability, the 
' energy flow and the full width at half maximum(FWIIM) of the surface solitons. We state that the refractive index 
difference affects the power distribution of the surface solitons in two media. In addition, we show that the different 
boundary values at the interface can lead to the different peak position of the surface solitons, but it can not influence 
the solitons stability with a certain degree of nonlocality. 

Here, we consider the simple (l-l-l)D case, a; = is the interface of the nonlinear medium(a linear refractive index 
n^) and a linear medium of refractive index ng. To describe the propagation of light beams along Z axis near the 
interface of the nonlinear medium, we use a nonlinear Schroinger equation for the dimensionless amplitude a of the 
light field coupled to the equation for normalized nonlinear induced change of the refractive index (j), 
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where X is is the transverse coordinate, = w'^/wq, wq is the beam width, Wm is the characteristic length of the 
nonhnear response and a stands for the degree of nonlocahty of the nonhnear response. /? and fco = 2'k/X are the 
wave numbers in the media and in vacmim. </> is given by (t>{X) — ^oe^/" — (l/o!^) /"^ G(X, ^)|a(^)p(i^, where 

G{X) = ae«/"(e^/" - e--^/")/2, for X > ^, and G{X) = ae-^/"(e«/" - e-«/")/2, for X < ^ Here, we can safely 
assume that the boundary condition at the interface(X = 0) is = 0o(0o > 0), where (/jq = 50 is the initial value, 
unless we indicate it. a and (j) vanish at the X — J> ±oo. 

We search for stationary soliton solutions of Eqs. (H]) and ([2]) numerically in the form a{X, Y, Z) = u{X, Y) exp{ibZ), 
where u is the real function and & is a real propagation constant of spatial solitons in the normalized system. 

1 w'^ 

-\/'iu~bu + ^{k'^nl- I3'^)u + (f>u = 0, (3a) 
1 i/;^ 

-Vlu-bu+^ {klnl -I3')u = 0. (3b) 

and 

a^Vi^-^+lup = 0. (4) 

To elucidate the linear stability of the solitons, we searched for perturbed solutions in the form a{X, Z) — [u{X) + 
p(X, Z) + iq{X, Z)] exp{ibZ) i. Hi [H, [H , where the real p{X, Z) and imaginary q{X^ Z) parts of the perturbation 
can grow with a complex rate 5 ~ 5r + iSi upon propagation. Linearization of Eq.(3)and (4) around a stationary 
solution yields the eigenvalue problem 

6p=-l^ + bq-H- %klnl - P')q, (5a) 

Sq^l-^-bp + cl>p+ ^(fc^ni - P^)p - uAcl>. (5b) 

which holds for -oo < X < 0. where Acj) = 2 G{X,^)u{^)p{^)d^. 
For < X < oo, the eigenvalue problem 

Sp='l^ + bq-f{kK-(3'H (6a) 

Sq=l^,^bp+fiklnl^P')q. (6b) 

We first consider that the influence of the difference of and uq on the surface solitons. For the zero field 
refractive index of the nonlinear medium(nL) is larger than the refractive index of the linear medium(no), at the same 
degree of nonlocality of the nonlinear response, from Fig. [TJ we can find that the solitons reside almost fully inside the 
nonlocal nolinear region and only weakly penetrate into the linear region when two media have a large refractive index 
difference, but the surface solitons have a significant part of their optical power residing in the linear medium when 
the boundary is between two media with a small refractive index difference which is comparable to the nonlinear index 
change. The results show that the refractive index difference affects the power distribution of the surface solitons in 
two media. Comparing Fig. [TJb) with (c), we can see that the profiles of surface solitons are alike when the refractive 
index difference between two media is small. However, when the refractive index difference between two media is 
big(Fig. [Ija) and (d)), the profiles of solitons are very different and solitons are no longer affected by nonlocality 
shown in Fig. [2l^d)(no — = 0.6). Of course, when the degree of nonlocality a is equal to zero, that is to say, the 
nonlinear media is local, the solitons are stable in the case of Fig. [Ijc), whereas the solitons are unstable in the case 
ofFig.lHb) 

The central finding in this Letter is the influence of the change of nonlocal degree on the solitons stability. In 
Fig. IHa), for > uq, with the degree of nonlocality becomes stronger, the solitons are more stable. When the 
degree of nonlocality exceeds a certain value, the solitons will be stable. The index difference influences the value. For 
n-L < no, only when the index difference is small, the solitons stability will be affected by the degree of nonlocality. 
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This is shown in Fig. [2jc). Fig. [UJb) and (d) depict the sohtons are very stable, propagating without distortion or 
deviations in their trajectories for a propagation distance of 15 diffraction lengths with 5% white noise. These results 
illustrate the fact that the nonlocal nonlinearity does action on the surface soltions. when the force exerted on the 
beam by the nonlocal nonlinearity is equal to the force exerted by the boundary at the interface, the solitons keep their 
straight line trajectories. Here, we only show that the cases — no = 0.6[Fig.[2jb)] and no — = 1 x 10^^[Fig.[2][d)]. 
Having demonstrated the influence of the degree of nonlocality on stability of the surface solitons, we proceed to study 
the energy flow P = \a\'^dX or FWHM of the surface solitons as a function of the degree of nonlocality a[Fig. 
As the degree of nonlocality increase, the energy flow monotonously increases. FWHM firstly increases with the 
increase of a, but it will decrease when the degree of nonlocality is strongly nonlocal. Importantly, the boundary 
value at the interface can also dramatically modify the properties of surface soltions. For example, it can affect the 
position Xjnax of the maximum value of \u\ [Fig. H^a)]. X^ax is located farther away from the interface when the 
boundary value 0o is smaller. In Fig. ^b) and (c), one can easily find this point by comparing the surface soltion 
at 00 = 20 with the surface soltion at (j)Q = 50. So, we can say that the force exerted on the surface solitons by the 
interface will increase when the boundary value increases. The force attracts the surface solitons to the interface. 
However, the boundary value at the interface can not infiuence the stability of solitons when a is a certain value. This 
can be explained by Fig. Sljd) in which the change of the perturbation growth rate Sr followed by (j)o is a straight line. 

To summary, the stability, energy fiow and FWHM of the surface wave solitons can be affected by the degree of 
nonlocality of the nonlinear medium. We find that the refractive index difference affects the power distribution of the 
surface solitons in two media. We state that the different boundary values at the interface can lead to the different 
peak position of the surface solitons, but it can not infiuence the solitons stability with a certain degree of nonlocality. 
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No. 06CXTD005). 
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FIG. 1: Profiles of surface solitons at a = 0, and a = 20 for (a)n_L — no = 0.6 , (b)n_t — no = 1 x 10 ^, (c)no 
and (d) no — ul = 0.6. 
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FIG. 2: The influence of the degree of nonlocaUty q on the perturbation growth rate 5r for n_L — no = 0.6 or til — no = 1 x 10~''(a) 
and no — nL = 1 X 10"'' (c). Propagation of the surface solitons launched at a; = with 5% noise at q = 20 for ul — no — 0.6(b) 
and no — ni = 1 X 10~''(d). 
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FIG. 3: Energy flow P and FWHM versus a for surface solitons for ul ~ no — 0.6(a) and no — = 1 x 10 ®(b). 
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FIG. 4: (a)The position Xmax of the maximum value of \u\ as a function of the boundary value (j!>o at the surface. Circles 
correspond to surface soltion at (f>o = 20 or <^o ~ 50 shown in (b) for ul — no — 0.6 and (c) for no — ul ~ I x W~^'. (d) The 
influence of <^o on 5^. a ~ 20 for all figures. 



